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^ . I. INTRODUCTION 

X) , 

' In recent years a isospin chemical potential has been introduced into the study for QCD phase structure [H, 0, 0] . 
[ It allows one to have a new dimension to study the rich phases of QCD theory. When compared to the baryon 
• chemical, this isospin chemical potential in principle has no fermion sign problem on a lattice simulation From 
recent lattice calculation and the investigation of chiral perturbation theory at finite isospin chemical potential, it has 
been indicated that there will be a pion condensation or pion superfluid phase 0i Q • At finite temperature and finite 
^ ' isospin density it is in essence a Bose-Einstein condensation (BEC) of charged pion in momentum space. This kind 

, phenomenon has been also studied by using NJL model [1,0,11 and ladder QCD [i]. 
CO ' In the low energy effective models, the linear sigma model is very simple and illustrative. When only considering 
the mesonic part of the model, there are only four scalar fields, the sigma field and the usually three pion fields which 
' display a 0(4) symmetry. The model is so called 0(4) linear sigma model. It has been usually used to study the 
chiral phase transition and it is also well suited for describing the physics of meson 121 1. This model has been 

. well studied at finite temperature within the Cornwall-Jackiw-Tomboulis (C JT) formalism jlO| . In our previous work, 
we have introduced the isospin chemical potential into the linear sigma model and discussed the BEC of pions and 
its relation to the chiral phase transition in the chiral limit within the same formalism [l3l |. In this paper, we will 
extend this work to discuss the BEC not only in chiral limit but also in non-chiral limit at the Hartree approximation 
Oh, and also at the large N approximation. Through this work we wish to give a clearer picture of BEC in linear sigma 
^ • model at finite temperature and finite isospin density. 

The organization of this paper is as follows. In section 2 we give a brief introduction of the 0(4) linear sigma model, 
then we introduce the isospin chemical potential. The CJT formalism will be briefly described based on Xcj}* theory. 
In section 3, the CJT method is used to derive the gap equations and thermodynamic functions, then we will discuss 
the BEC in chiral limit and non-chiral limit at Hartree approximation. In section 4, the BEC is discussed at large N 
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' approximation. The last section is the summary. 

II. THE LINEAR SIGMA MODEL AT FINITE ISOSPIN CHEMICAL POTENTIAL AND THE CJT 

FORMALISM 

We start our discussion from the Lagrangian of the linear sigma model with only the mesonic part presented, 

- ^(^^ + 7r2)^+ea, (1) 

where a and 7? are the sigma field and the three pion fields (tti, 7^2,1^5) respectively, ea is the explicit chiral symmetry 
breaking term, where e = f-wi^t and /tt — 93MeV is the pion decay constant. At tree level and zero temperature 
the parameters of the Lagrangian are fixed in the way that the masses agree with the observed value of pion mass 
niTr = 138Mey and the most commonly accepted value for sigma mass nia = GOOMeV. Then the coupling constant 
A and negative mass parameter of the model are chosen to be A = 3(m^ — mj)//^ and — = (m^ — 3m^)/2 > 0. 
When e = 0, the chiral symmetry is spontaneously broken, and pion is the Goldstone boson which mass is TOtt = at 
zero temperature. 
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The isospin chemical potential can be introduced by different means. In Ref. '14'|, the chemical potential is intro- 
duced according to the covariant way fixed by the gauge invariance. In our previous work |.13i] , the chemical potential 
is introduced through the conserved isospin charge. Both give the identical results. If we redefine the pion fields as 

TT- = -^{tti +iTT2), 7r+ = -^(tti - i7r2), 7ro = 7r3. (2) 



the Lagrangian with the isospin chemical potential /i included can be written as 

1 
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where the potential is 



V{a, n) = ^(a^ + n') + + n') ea, (4) 

and (Dqtt)^ = {do ± ifi)n±, {DiTr)"^ = din± for i = 1, 2, 3. 

By shifting the sigma field a.s a ^ a + <j), where (j> is the expectation value of the sigma field and it is also the order 
parameter of the chiral phase transition, the classical potential takes the form 

while the interaction Lagrangian which describes the vertices of the shifted theory is given by 

Lint = a TT (T TT (pan (per . (6) 

12 24 24 6^ 6^ ^ ' 

By a generating functional theory, one can obtain the effective potential based the above Lagrangian. At finite 
temperature the effective potential is identical to the thermodynamic potential which is very important in discussing 
the thermodynamic properties of the system. The usually effective potential V^((/>) depends on 0, a possible expectation 
value of the quantum field $. A generalized effective potential for composite operators has been introduced by 
Cornwall, Jackiw and Tomboulis(CJT) TE\. According to this formalism, the effective potential V{(t), G) depends not 
only on 0, but also on G{x,y), a possible expectation value of the time-ordered product T^{x)^{y), which is also the 
full propagator of the field. Physical solutions demand minimization of the effective potential with respect to both (jj 
and G, which means 

dcP ' dG ^ ' 

The derivations to (j) and G are functional. This formalism was originally written at zero temperature. Then it was 
extended to finite temperature by Amelino-Camelia and Pi for investigations of the effective potential of the 
theory 

In our following discussions of linear sigma model at finite temperature, we use the imaginary time formalism, 
which is also known as the Matsubara formalism [17]. This means 

-'^Y.j |^/(^-". « - / («) 

where (3 is the inverse temperature, (3 = 1/T; the integration over the time component kg has been replaced by a 
summation over discrete frequencies. For boson there are tUn = 2TmT and n = 0, ±1, ±2, • • • . For the sake of simplicity 
in what follows, a shorthand notation is used to denote the integration and the summation. 

For At/)** theory, according to the CJT formalism [l^, the effective potential can be written as 

V{(l>, G) = U{4>) + i ^ lnG-i(0; fc) + i ^[^"'('/'; k)-l] + ^2(0, G), (9) 

where D and G are the bare and full propagators of the shifted A^* theory respectively. The last term V2((/), G) 
represents the sum of all two and higher-order loop two-particle irreducible graphs of the theory with vertices given 
by the interaction Lagrangian and propagators set equal to G(0;A;). The diagram contribution to V2((/', G) for Xcj)^ 
theory are shown in figure [T] </> and G will be self-consistently determined by equation ([7]). In the case of At/)'* theory, 
the first diagram of figure [IJa) is the leading order in V2(0, G) in both the loop expansion and the 1/iV expansion. In 
our later discussion of linear sigma model at Hartree and large N approximation, one needs to take into account the 
"00" type diagram only. 





(a) 






(b) 



FIG. 1; Two-particle irreducible graphs which contribute to the efTective potential of a Xcp"^ theory in the CJT formalism up to 
(a) two loop level and (b) three loop level. The solid line represents the full propagator. 



III. BEC AT HARTREE APPROXIMATION 



From the Lagrangian ([3]), after shifting the sigma field we can write down the tree level inverse propagators of a, 
ttq and t:± respectively as 



DZ 



k2 



D 1 = (cj„ + i^)2 + k2+m2 + -(/)2. 
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(10) 

(11) 
(12) 



The equation ()12p represents the inverse propagator of and 7r_. As the summation in equation (jSj) is symmetric 
over n from — oo to +cx3, uon + and uJn — i/i are equivalent in describing the propagators of 7r+ and 7r_. According to 
the CJT formalism, the effective potential at finite temperature and finite isospin chemical potential can be written 
as, 

y (0, G) = {/(0) + lj\nG-' + \j^ [D-'G, - 1] 

+ \j\nG^' + \j^[D^^G^-l]+j\nG-'+j^[D-'G-l]+V2{ct^,G), (13) 

where G^ , Gq and G are the full propagators of a, ttq and 7r± respectively. They are determined by the stationary 
condition ([7]). V2 (</!), G) represents the infinite sum of the two particle irreducible vacuum graphs. However, at the 
Hartree approximation we need only to calculate the "00" (or "double bubble" ) diagrams and treat each loop line as 
the full propagator [l^, [l^. Therefore, V2 can be written as 



^2(0, G) 
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For the full propagators we could take the following ansatz. 



G-i = 



k^ + 



G^^ = ujI+\^ 



a 1 

Ml 



G-i = (cj„ + iM)'+k2 + A/2, 



(14) 



(15) 
(16) 
(17) 



where the effective masses Mcr,Mo and M have been introduced for ct, ttq and Tr± respectively. As in the Hartree 
approximation, only the tadpole diagrams contribute to the self-energies. The effective masses are independent of 
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momentum. From the stationary condition ([7]), we obtain a set of effective mass gap equations, 
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Accordingly the effective potential can be written as 

V{4>.M) = [/(<^) + i / lnG;i 



^ / Go. 



(18) 
(19) 
(20) 



[Ml 



m 



)Go 



InGn 



{M^ -m" -^4>'')G + V2{(j), M). 
By minimizing the potential with respect to the order parameter 0, we obtain one more equation, 
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(21) 



(22) 



From equations (fT5|) — (|20p and (P^ . the effective masses and order parameter can be solved self-consistently at given 
temperature and chemical potential. There are four types of integrals in above equations. After performing the 
Matsubara frequency sums they are 

Q(M) = / ln[{uj„ + ifif + + M2] 



MM) + Qf}iM) 



1 



1 



F{M) 



[3 J (27r)3 2 P J (27r)3 
1 

'p {Ujn + i/i)2 + k2 + Af 2 

F„{M)+Fp{M) 



ln(l 



(23) 



(fk 1 



(27r)3 2w 7 (27r)3 2^ 



d^k 1 



1 



1 



(24) 



Q{S) = / ln(cj2 +k2 + 52) 



)o(5) + g/3(5) 



d^k 



F(5) = 



1 

(27r)3 2 ' /? 
1 

^ a;2 + k2 + ^2 

i^0(^) + i^/3(^) 
d^k 1 



(27r)3 2wk 



d3k 

(2^ 



d3k 1 



ln(l - e-'^"-). 



(25) 



(27r)3 Wk e'^'^'' - 1 ' 



(26) 



where lo — V^ + AP , Wk — v^+S^ and 5 = M^, Mq. Each integral is divided into two parts: the zero-temperature 
part which is divergent and the finite- temperature part which is finite. The evaluation of the integral requires 
renormalization. There are some discussions concerning the renormalization on the CJT formalism in different mod- 
els [l^[l^,[l^. The investigations about the renormalization of the O(A^) linear sigma model can be referred to [Tll.[20j. 
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In our discussion, we only kee p th e finite temperature parts {Qp and F/j). The divergent parts of the integrals are 
neglected as is done in O, [11,111. 

In the discussion of the thermodynamic system, the effective potential V is equivalent to the thermodynamic 
potential Q. Thus we have 



n 



1 



1 
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A 



6 



(27) 



where 



= ^mAU)f + ^[F,iM)f + ^[Fp{Mo)f 



^Fp{M„)Fp{M) + ^Fp{M,)Fp{M^) 



A 



Ffi{M)F0(Mo). 



(28) 



According to the relation 



(29) 



and equation ([20|) . we can get the net charge density 



(2^ 



1 



g/3(w+/j) 



1 



(30) 



This expression of density seems very similar to that of ideal gas, but actually they are different. Because here the 
effective mass M in a; is a function of temperature and will be determined self-consistently by the gap equations. 

Now we are in a position to discuss BEC. If we lower the temperature of the system from a high temperature, it is 
known that BEC will occur at the place where ji — M \s reached. When BEC occurs, equation ([50)1 should be written 
as 



(27r)3 



1 



1 



(31) 



where po represents the charge density of the zero-momentum state [ij 



A. Chiral Limit e = 



In the chiral limit {e = 0), the corresponding coupling constant and the negative mass parameter are given by 
A « 125 and — « 1.8 x lO^MeV^. From equation (|30p . when p is fixed, by solving the gap equations ^5)1 — ([^0]) 
and equation ([22]), we find both p and M are functions of T. It can be plotted in figure [2l We can see that with 
T decreasing, M decreases first and then increases, while p keeps increasing quickly and approaches M. At certain 
temperature, p catches up with M and BEC happens. From equation (31), we know the critical temperature is 
determined implicitly by the equation 

p^p*{l3,,p^M). (32) 

When T < Tc, the system goes into BEC phase. The equation (pij) will be solved together with the gap equations 
with p fixed. We find p and M are still functions of T and p{T) = M{T). They both decrease with temperature 
decreasing, which is indicated in figured 

For different fixed total density p, there will be different values of critical temperature Tc and chemical potential 
p(Tc), so one can study the p — T phase diagram of BEC. It is shown in figure [31 In the chiral limit, at T = 0, the 
pion mass M = 0, so the curve starts from the point T ~ and p = 0. There is a disjunction at T « lASMeV. 
This is due to the first order chiral phase transition in the chiral limit at Hartree approximation. At low temperature 
(/) ^ 0; at certain high temperature = 0. From (/) ^ to = is discontinuous, which is reflected in the BEC 
phase diagram as a disjunction from low temperature to high temperature. Thorough discussion about the relation 
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FIG. 2: ^ and M as functions of T at the fixed total charge density p = 0.06/m ^ . BEG happens at Tc — SOMeF. 

of the BEC and the chiral phase transition in the chiral Umit can be referred to our previous work [l^ . In the chiral 
symmetry broken state (</> 7^ 0), the pion mass is non-zero {ji = M ^ 0) which means the Nambu-Goldstone(NG) 
theorem is not observed at Hartree approximation as also indicated in the previous literatures [10, 21, 22]. Recent 
discussions of this direction can be referred to [H,!!^!. In this paper, we mainly discuss BEC at the usually Hartree 
approximation and large N approximation. In our later discussion, we will see that at large N approximation the NG 
theorem will be preserved. In figure [H we can see clear that the whole phase plane has been divided into the BEC 
phase and the normal phase. 
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FIG. 3: The phase diagram of /i versus T for BEC in the chiral limit at Hartree approximation. 



B. Non-Chiral Limit £ / 

In the non-chiral limit (e 7^ 0), the corresponding coupling constant, the negative mass parameter and the explicit 
symmetry breaking term are given by A « 1.18 x 10^, — « 1.51 x IQ^MeV"^ and e « 1.77 x IQ^MeV^. When 
equation ([50]) is solved together with the gap equations (fTS]) — (|20l) and equation at fixed p, the critical temperature 
Tc of BEC will be determined aX ^ — M which is shown in figure |4l The procedure is similar to that of the chiral 
limit. For different fixed p, /i varies with Tc, so the /i — T phase diagram of BEC in the non-chiral limit can be 
plotted as shown in figure [5] At zero temperature, the vacuum mass of pion is IZ^MeV, so the curve starts from the 
value of /ic — JTT-TT — li^MeV . The chiral phase transition in the non-chiral limit is a smooth crossover. There is no 
discontinuity in the BEC phase diagram. 
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FIG. 4: fj, and M as functions of T at the fixed total charge density p = 0.3 fm . BEG happens at Tc — lS2MeV. 
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FIG. 5: The phase diagram of jj. versus T for BEG in the non-chiral limit at Hartree approximation. 



IV. BEC AT LARGE iV APPROXIMATION 

The generalized version of the meson sector of the linear sigma model is called 0{N) model and is based on a set 
of N real scalar fields. The 0{N) model Lagrangian can be written as 

C = - im2*2 „ _La*2 + (33) 

where $ can be identified as $ = (a, 7ri,7r2, • • • ,7rAr_i). The last term is the symmetry breaking term in order to 
generate the observed masses of the pions. If = 4, it becomes the 0(4) linear sigma which has been discussed 
above. 

For the 0{N) linear sigma model, after introducing the chemical potential and shifting the sigma field the Lagrangian 
can be written as 

o\ 1 1 

i 

_ i(^2 ^^2) y 2 _ 1(^2 ^^2)^2 _ 1^2^2 „ ^ 

i 

where |7rp = tt+tt-, the sum over i or j is from 3 to iV — 1 and i ^ j- The chemical potential here is associated with 
the conserved charge of an 0(2) symmetry [25'|. Then we can write down the tree level inverse propagators of tr, tt^ 
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and 7r± respectively as 



,2 I 2A 2 



N 
2A 

3iV^ 



{uJn + illf 



2A 



(35) 
(36) 
(37) 



From the CJT formalism, the thermodynamic potential of the 0{N) linear sigma model can be written as 



+ 



f)(0,G) = imV + ^0 

iV-3 



1 



1 



- / lnG-i + - / \D-^G,-l 
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^ 3N 
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(38) 
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(39) 



By minimizing the thermodynamic potential with the full propagators, we obtain the following set of effective mass 
gap equations 



Ml 
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Gi, 



(40) 
(41) 
(42) 



P 



where Mi stands for the effective mass of tt^ (i = 3, • • • , TV — 1). By minimizing the thermodynamic potential with 
the order parameter 0, we have one more equation 



2A 



2A 



3iV 



iV 



2(7V-3)A 
3iV 



e = 0. 



(43) 



In the large iV approximation, which means that we ignore the terms of 0{1/N), the equations PO]) — (|43p reduce 



to 



M^=™2 + ^<^^ + ^F,(M,), 
= m2 + |A<^^ + ^F,(M,), 

M^ = + ^ '^p Y 

3N^ 3 '^^ ' 



^2 + (M, 
3A^^ 3 '^^ ■ 



e = 0. 



(44) 
(45) 
(46) 
(47) 



The terms quadratic in (}> are not of 0{1/N) but of 0(1) since (j) depends on N as, (fP' — 3Nm^/2X. For the integration 
only the finite temperature part is preserved as done in Hartree approximation. The effective mass equations (|45p 
and (I46p are identical which shows the effective mass of 7r± is equal to that of tt^. 
From the thermodynamic potential we can also derive the net charge density as 



(2^ 



1 



1 



(48) 



where uj = Vk^ + M^. In the following discussion we will set iV = 4 for our numerical evaluation. 
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A. Chiral Limit e = 

In the chiral limit, by the equation (|46)) . the equation ([47)1 could be written as 

(j)M^ = 0. (49) 
When 4> ^ 0, the effective mass of pion will be zero. Therefore the pions are massless Goldstone bosons in the chiral 

300 
250 
200 

I "° 
100 

50 



50 100 150 200 250 300 

T (MeV) 

FIG. 6: fj, and M as functions of T at the fixed total charge density p = BEG happens at — 190MeV. 

symmetry broken phase which means the NG theorem is observed in the large N approximation. In chiral symmetry 
broken phase BEC will happen a,t fi = M = 0. At certain high temperature that the chiral symmetry restored and 
= 0, pion becomes massive because of the thermal contribution to the effective mass. In this case the equation p6)) 
and (|48p will be solved together at fixed density to determine the critical temperature of BEC as shown in figure [S] 
For different density by determining the critical temperature and chemical potential, we can plot The — T phase 
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FIG. 7: The phase diagram of fi versus T for BEG in the chiral limit at large A'^ approximation. 

diagram of BEC as show in figure [7l At low temperature when chiral symmetry is not restored, the BEC will happen 
at /i(r) = M(T) = which reflects the requirements of the NG theorem; At high temperature when chiral symmetry 
is restored, the BEC happens when fi equal to the nonzero effective mass of pion. 

B. Non-Chiral Limit £ / 

In the non-chiral limit, the equation ([46]) . (|47p and ([48]) will be solved together at fixed density to find the critical 
temperature of BEC. The figure [5] shows the critical temperature at the fixed density is determined just at the time 
/Lt(T) = M{T). By the same procedure as that in the chiral limit, the fi — T phase diagram could be plotted as in figure 
ini At zero temperature the vacuum mass of pion is 138MeV, and the BEC happens at critical chemical potential 
He = I'iSMeV. When temperature increases the critical chemical potential also increases. The BEC phase is in the 
upper plane of the phase diagram. 
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FIG. 8: fx and M as functions of T at the fixed total charge density p = 0.1/m~^. BEG happens at Tc = 125MeV. 

600 
500 

400 
I 300 
200 



100 



50 100 150 200 250 300 350 400 
T (MeV) 

FIG. 9: The phase diagram of fj, versus T for BEG in the non-chiral limit at large N approximation. 



V. SUMMARY 

By the CJT formalism we have derived the temperature and density dependent effective potential based on the 
linear sigma model. The BEC is investigated at the Hartrcc approximation and the large approximation. The 
critical temperature of BEC is determined by lowering the temperature of the system at the fixed density to find the 
critical point at which = M(T). The fi — T phase diagram of BEC has been plotted in different situations. The 
main features of BEC at different approximations can be summarized as below: 

(1) At Hartree approximation and in chiral limit, for the /z — T phase diagram of BEC, the phase separate line 
is discontinuous from low temperature to high temperature. At low temperature phase {(f) ^ 0), BEC happens at 
^ = M > 0, which shows the NG theorem is not observed. 

(2) At Hartree approximation and in non- chiral limit, the critical chemical potential of BEC at zero temperature 
is exactly the vacuum mass of pion (138MeF). With temperature increasing the critical chemical potential increases 
continuously. 

(3) At large N approximation and in chiral limit, at low temperature phase before chiral symmetry restored 

((/) 7^ 0), BEC happens exactly at /i(r) = M{T) = 0, which shows the NG theorem is observed. Above certain 
temperature when chiral symmetry is restored, the critical chemical potential of BEC becomes non-zero and increases 
with temperature increasing. 

(4) At large N approximation and in non-chiral limit, the critical chemical potential of BEC is 138MeV at zero 
temperature then increases with temperature increasing. 
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